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We investigate the impact of the matter and glue dynamics on baryon number fluctuations and 
the kurtosis of baryon number distribution. This is done within the framework of QCD-improved 
low energy effective models. In particular we include the momentum scale dependence of the quark- 
meson scattering and the non-trivial dispersions of both, quarks and mesons. On the gluonic side 
we take into account the backreaction of the matter sector on the glue dynamics. It is shown that 
the above fluctuations lead to a more rapid change of the baryon number fluctuations as well as the 
kurtosis of with the chiral crossover. 

We also study the signatures of quark confinement in low energy QCD. It is shown that contrary to 
the common picture the effective thermal distribution in the presence of confining glue backgrounds 
does not tend towards the colourless baryonic one. Instead, the dominance of colourless hadronic 
states is obtained in a subtle interplay of quark and glue contributions to the canonical potential. 
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I. INTRODUCTION 

Studies of the QCD phase structure and the QCD ther¬ 
modynamics have attracted lots of attentions in recent 
years. QCD matter composed of deconfined, coloured 
quarks and gluons, i.e. the Quark-Gluon Plasma, pro¬ 
duced before the quark-hadron transition at ~ 10 /iS in 
the evolution of early Universe [T] , is believed to be recre¬ 
ated in heavy-ion collisions at Relativistic Heavy-Ion Col¬ 
liders such as (RHIC) [2ll3] and LHC [4j. 

The chiral and confinement-deconfinement crossovers 
for QCD with physical quark masses at vanishing density 
may turn into first order transitions at high density. One 
of the key challenges concerning the phase structure of 
QCD is to get hold of the existence and location of the 
corresponding critical end point (CEP), [5]. The Beam 
Energy Scan (BES) program at RHIC is directly aimed 
at this task, where the beam energy or collision centrality 
dependence of the fluctuations of conserved charges, such 
as moments of net-proton multiplicity distributions, is 
employed to locate the QCD critical point. 

In turn, the experimental programme should be accom¬ 
panied by reliable theoretical predictions for the above 
observables and their relation to the CEP are highly de¬ 
manded. A particularly promising direction is the com¬ 
putation of the fluctuation of conserved charges at van¬ 
ishing and finite density, see e.g. [6H9]. Moreover, it 
has been suggested that the values of these observables 
at small density already give access to the existence and 
location of the CEP, see e.g. m- 

In the present work we investigate the QCD phase 
structure and thermodynamics as well as baryon num¬ 
ber fluctuations within the quark-meson (QM) model 
and the QCD-enhanced Polyakov-quark-meson (PQM) 
model. The thermal and quantum dynamics of the mod¬ 
els is accessed with the functional renormalisation group 
(ERG), for QCD-related reviews see [TTHI^ . Recent de¬ 
velopments in ERG applications to QCD, and in partic¬ 
ular the embedding of low energy effective models in first 


principle QCD can be found in [161 USHlS], for a recent 
survey see [27j . Baryon number fluctuations and general 
higher moments have been studied with the ERG in |5S]- 
15^ . see for corresponding mean field results and 

interesting algorithmic developments. 

Here we build on the progress made in QCD-embedded 
low energy effective model, see [IllIIHil, to improve 
on the existing fluctuations studies within an advanced 
ERG study of these observables: the momentum scale de¬ 
pendence of the quark-meson scattering is taken into ac¬ 
count as well as the non-trivial dispersions of both quarks 
and mesons. Furthermore the backreaction of the mat¬ 
ter sector on the glue dynamics is included, leading to 
a temperature-dependence modification of the Polyakov 
loop potential. 

The paper is organized as follows. In Section[Tl]the em¬ 
bedding of low energy effective models in QGD within the 
ERG framework is recalled. In Section [TTT| the QM model 
is discussed. The significance of the additional mat¬ 
ter fluctuations considered in the present work is most 
cleanly seen without the Polyakov loop potential. We 
derive the flow equations for the couplings, and present 
numerical results on the phase structure, thermodynam¬ 
ics and baryon number fluctuations. In Section |IV| we 
discuss the QCD-enhanced PQM model. This includes 
a discription of the QCD-enhancement of the Polyakov 
loop potential, as well as an discussion of the intricacies 
of the relation between color confinement and hadronic 
properties. At the end of this Section our final results 
are presented for the order parameter, the thermody¬ 
namics and the higher moments including a 2-1-1 flavour 
estimate for the kurtosis in comparison to lattice data. 
Some useful formulae, such as threshold functions, flow 
equations, and further discussions can be found in the 
Appendix. Notably, in Appendix we discuss the fre¬ 
quency and chemical potential dependence of fermionic 
couplings. This is done in view of the silver blaze prop¬ 
erty, and the thermal decay properties specifically rele¬ 
vant for the thermodynamics at low temperatures. 
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FIG. 1. Flow equation for the effective action or grand po¬ 
tential F or n in full QCD. The first two loops with gluons 
and ghosts give the glue contribution flgiue, the quark and 
mesonic loops give the matter contributions flmatt- 
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FIG. 2. Flow equation for the effective action or grand po¬ 
tential r or in the low energy regime. The gluon and quark 
loops in Fig. are suppressed and the quark and mesonic 
loops give the matter contributions fimatt. 


II. FROM QCD TO LOW ENERGY EFFECTIVE 
MODELS 


It has been argued, see e.g. m, that low energy ef¬ 
fective models can be derived systematically from QCD 
allowing for a precise determination of the fundamen¬ 
tal model parameters. Within the functional renor¬ 
malisation group framework this can be discussed al¬ 
ready on the level of the flow equation for the scale- 
dependent effective action, with the super-field 

$ = c, c, q, q, </>,...) with (p = (cr, 7?) and possible fur- 
the effective hadronic fields. The flow equation reads 
schematically 


dtTkm = ^TiG^^mdtRt, t = \n{k/A), (1) 


where 


G<s>^,s>j [*&] 


<5Tfc[4>] 

, ( 5$2 


Rt, 


( 2 ) 


is the full field-dependent propagator, k is the infrared 
cutoff scale, and A is some reference scale. The flow equa¬ 
tion 0 is depicted in Fig. for more details and QCD 
results, see e.g. uninsMTi. All loops in the flow carry 
only momenta q^ < with the infrared cut-off scale k, 
and 0 implements a successive integration of momen¬ 
tum modes. The first two terms in Fig. constitute the 
contributions of the glue system, while the third term 
entails the quark fluctuations and the last loop stands 
for the fluctuations of the effective hadronic degrees of 
freedom. Hence, the effective action can be written as 


rfc[$] + 


natt,fc 


— ^q,k + r 




( 3 ) 


where includes the integrated flow of gluon and 

ghost loop, that of the quark loop, and 

that of the hadronic degrees of freedom. Note, that the 
split in quark and hadronic contributions does not re¬ 
flect an effective theory setup. Within the framework 
of dynamical hadronisation, introduced in [HEHIS], it 
has turned out to be a very effective and powerful pa- 
rameterisation of matter fluctuations in ab initio QCD 
in terms of genuine quark scatterings and resonant mo¬ 
mentum channels with hadronic quantum numbers, for 
applications to QCD see e.g. pSl 1^ . 


Moreover, it facilitates the embedding of low energy ef¬ 
fective theories of QCD such as NJL-type and QM-type 
models with and without confining Polyakov loop poten¬ 
tial in first principle QCD: For the sake of definiteness 
we restrict ourselves to Landau gauge QCD. There, the 
physical mass gap of QCD is reflected in a mass gap in 
the gluon propagator of about 1 GeV. Accordingly, glue 
fluctuations decouple below this scale. The ghost has 
no mass gap but only couples to QCD matter via the 
gluon. Hence, it effectively shares the gluon decoupling. 
In summary, the fluctuations of the glue sector decouple 
from the matter sector in the infrared. Still, the effective 
potential of the glue sector, is required for 

solving the quantum equations of motion for low ener¬ 
gies. There, it is relevant for determining the physical 
glue background. In summary, integrating the QCD flow 
down to scales k = A with A < 1 GeV leaves us with 
an effective matter theory with glue background. The 
flow equation for this QCD-embedded low energy effec¬ 
tive theories is given by the two matter loops in Fig. 
see Fig.[^ The total effective action is given by 0 where 
Tgiuo.fc is an external input. Since it does not enter the 
flow equation of the matter part, we only need it at fc = 0. 

In a flow equation approach to QCD, Hgi^e comprises 
the contributions from the gluonic and ghost loop in 
Fig.[^ while comprises that from the quark 

and mesonic loop in Fig.[^ The representation in terms 
of a flow equation facilitates the discussion of the repsec- 
tive scales. Dropping the grand potential reduces 
to the pure glue potential with a transition temperature 
Tgiue ~ 200 MeV, see (THl HU |52] . The glue potential (in 
the current approximation) features a first order phase 
transition, and below the Polyakov loop expectation 
vanishes. The possible minor center-breaking corrections 
of the glue propagation due to the vacuum polarisation 
are neglected. We shall see, that for the present argu¬ 
ment we can safely drop these contributions as they only 
increase the expectation value of the Polyakov loop. 

We conclude that the non-vanishing expectation val¬ 
ues of L and L in QCD are dominantly triggered by the 
quark loop in The influence of the non-trivial Aq- 

or L, L-background will be discussed at the end of this 
work. For the time being we concentrate on the impact 
of improved matter fluctuations on the higher moments. 
This concerns in particular the influence of the deformed 
dispersion that relates to non-trivial wave function renor¬ 
malisations of quarks and mesons. 
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III. QUARK-MESON MODEL 


In the spirit of the arguments of the last section we first 
concentrate on the low energy matter quantum, thermal 
and density fluctuations for two-flavour QCD. These fluc¬ 
tuations are well-captured within a quark meson model 
at scales well below 1 GeV, see the reviews and 

references therein. Moreover, full QCD calculations in¬ 
dicate that the width of the transition region from the 
quark-gluon regime to the hadronic one is small. This 
guarantees that the effective model is applicable and re¬ 
liable at low energy, see EnUMS]. The scale-dependent 
effective action of the quark-meson model is given by 




+ hkq (r^a + ■ lij q + t4(p) - C(t| -|- ■ 


(4) 


with dxo J d^x, and the dots indicate higher 

order terms in derivatives and fields. Accordingly, Q 
has to be understood in the spirit of both, a derivative 
expansion at low energies as well as a vertex expansion 
in terms of the fluctuation physics. The validity of the 
latter is important for having a good grip on the fluctu¬ 
ation physics at low energies, and is crucial for a good 
access to higher moments. Note that the validity of the 
expansion in terms of the fluctuation physics is based on 
the relevance ordering of fluctuating modes that propa¬ 
gate off-shell in the loops. This should not be confused 
with the relevance of higher asymptotic states in the ef¬ 
fective action. This important distinction is at the root of 
the quantitative agreement for the QCD thermodynam¬ 
ics computed in Polyakov quark-meson models including 
the trace anomaly, see |22j . In turn, in approaches based 
on an expansion in asymptotic states a large number of 
the latter is necessary to have quantitative access to the 
full fluctuation physics. 

The meson field (j) = {a, tt) introduced in Q is in 
the 0(4) representation, with p = (/>^/2. Here k is 
the infrared (IR) cutoff scale of FRG, see Sec III A for 
more details; p, is the quark chemical potential. T 
are the SU{Nf) generators with Tr(T*T-^) = and 
]^0 _ 


= '^NfxNf The field-dependent effective po- 

V — 

tential 14(p) is 0(4) invariant, and the chiral symmetry 
is explicitly broken by the linear term —ca. Thus, the 
mass of the Pions tt is proportional to the linear break¬ 
ing parameter c. 

The quark and meson wave function renormalisations, 
i.e. Zq^k and are scale-, momentum-, and field- 

dependent. Moreover, at finite temperature Z^s split 
into Zll and Z-^, corresponding to those longitudinal and 
transversal to the heat bath, respectively. The transveral 
ones dominate the fluctuations physics. It has also been 
shown, that the approximation with scale-dependent Z’s 
already captures the effects of the non-trivial momentum 


Label 

Truncations 

LPA 


= 0, 

dth 

= 

0 

LPA' 


4 0, 

dth 

= 

0 

LPA -b hk 

^t^tpjq 

= 0, 

dth 

4 

0 

full 

dtZ,pfq 

4 0, 

dth 

4 

0 


TABLE I. Four different truncations and their labels used in 
this work. 


and frequency dependence for the propagators quantita¬ 
tively, see [24] . Hence, as a good approximation we only 
keep the scale-dependence for the wave function renor¬ 
malizations, and assume Z^^ = Z-^ throughout this work. 
Quark and meson fields interact with each other via a 
scale-dependent Yukawa coupling hk in Eq. Q. For the 
fermionic parameters Zq and h a consistent treatment of 
the frequency- and p-dependence is crucial for the cor¬ 
rect low temperature physics. This is discussed in detail 
in Appendix]^ Higher order quark-mesonic interactions 
can be included through the dependence of hk(p) on the 
meson field [23] . which will be discussed elsewhere. 

We investigate the thermodynamics of the QM model 
within various truncations. This sheds light into the roles 
played by the different fluctuations. Here, we denote the 
truncation with a full effective potential and both run¬ 
ning ZkS and hk, as ’full’. If the flows of Z’s and h are 
turned off, i.e. 

dtZ,p/q = 0, dth = 0, (5) 

with t in Q and the reference scale A being the ultra¬ 
violet (UV) cutoff, the effective action in Eq. Q is that 
of the local potential approximation (LPA). Apart from 
LPA we also investigate the LPA' approximation, i.e. the 
LPA with running Z’s, but with a constant h. Here we 
choose a constant renormalised Yukawa coupling, dth = 0 
with h = hjZ^^'^Zq, for more details see the next section. 
A vanishing flow of h provides a good approximation to 
the small scale-dependence of h in full QGD, see pEll^ . 
In turn, dth = 0 does not work well due to the large 
anomalous dimension of the mesons. Another approx¬ 
imation is characterized by a running hk but constant 
Z’s, which is denoted as LPA + hk in our work. In Ta¬ 
ble H] we summarise the truncations studied in this work. 


A. Flow equations 

As discussed in detail in Section [IT] the flow equation 
for the quark meson model follows from that of first prin¬ 
ciple QCD for low cutoff scales k < IGeV, where the 
gluonic fluctuations decouple. Then the matter loops 
carry the low energy fluctuations in the presence of a 
glue background. The glue sector only gives rise to a 
background potential for Aq or the Polyakov loop L,L 
respectively. Dropping (the flow of) the background po¬ 
tential and restricting ourselves for the time being to triv¬ 
ial backgrounds Aq = 0 or L, L = 1 we arrive at the flow 
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equation for the quark meson model, see Fig. 

= ^TTG^4^]dtRt - TTGgg[<i>]dtRl, (6) 

where the super field now only includes the fluctuating 
quark and meson fields, $ = (g,g, </>), and the regula¬ 
tor = {Rl , Rf) suppresses infrared fluctuations of 
quark and meson fields respectively. The traces in Eq. (1 m 
sum over momenta and internal quark (Dirac and color) 
and meson (flavour) indices. The full, field-dependent 
propagator, (§ reduces to the coupled one of quarks and 
mesons, to wit, 

G^rj,/qqm = f ^ J C^) 

\ k / (p(j)jqq 

for the diagonal parts that show up in (|^. 

B. Flow equations for the effective potential 

In this work we use 3d flat regulators, m iZ], for 
quarks and mesons, see also Appendix The flow of 
the effective potential 14 (p) is obtained by substituting 
Eqs. (|^,([^ in ^ and evaluating the flow at constant 
mesonic fields, 

dtVkip) = ^[{Nf - 


for rjq. The latter also contains an evaluation of the p and 
T dependence important for the silver blaze property of 
QCD and the thermodynamics at low temperatures. 

In general it is more convenient to work with renor¬ 
malized fields and renormalization group (RG)-invariant 
quantities. We denote them by symbols with bar, as the 
renormalized masses shown in Eqs. §. The relations 
between renormalised quantities and original ones read 

= ^ (11) 

taking the meson field and Yukawa coupling as exam¬ 
ples. Thus we have p = Z^^kP and the effective potential 
14(p) = Vkip)- Note that the chiral symmetry breaking 
term —ccr in Eq. ([^ is linear in the fields. Hence it nei¬ 
ther contributes to the right hand side of any flow, nor 
does it flow. Considering Ck = this leads to 

— '^Pip^kC-k ■ (1^) 

In the present work, we solve the flow equation for the 
effective potential, Eq. ([^, within a Taylor expansion 
about a fixed unrenormalised field value k, to be con¬ 
trasted to one about the scaling minimum po.fej see e.g. 
|48) . It has been shown in [53] that such a Taylor ex¬ 
pansion about a fixed k has the most rapid convergence. 
Expanded in renormalised fields, the effective potential 
then reads 


-I- ^(-8,4) 2 


(7^k 5 


P<t>,k\T) 


-m,Nflf'^\rhl^,Pq,k-,T,p)\ , ( 8 ) 

where e.g. [5011531, are the threshold functions, 

see Appendix]^ The renormalised dimensionless meson 
and quark masses are given by 


^2 


'■Ti-.fc 


vL{p) 

k'^Z, 


4>,k 


^2 


Vlip) + 2pVl'ip) 


^a,k 




4>,k 


^q.k 


‘^k^Zlk 


and the anomalous dimensions are defined by 

^t^<p,k _ ^t^q,k 


V4>.,k — 


■^(j>,k 


Vq,k — 


-'q.k 


(9) 


( 10 ) 


Here we consider frequency and spatial momentum- 
independent anomalous dimension, the anomalous di¬ 
mensions are evaluated at low frequencies and spatial 
momenta. The related depencence if covered by the k- 
dependence. The precise definition of the anomalous di¬ 
mensions in ( [To] ) and a detailed discussion of this approxi¬ 
mation is found in Appendix [Bjfor and in Appendix [Cj 


N 


Vk{p) = Z, 

n\ 

n—Q 


with \n^k = ^n,k/Z^^f. and Kk = Z^n with 
^t^k P(p,k^k ■ 


(13) 


(14) 


In accordance with the convergence discussion in |23| we 
already find a good convergence for N = 5. Eor example, 
the difference for /^ = (t between N = 5 and A^ = 6, 7 is 
less than 1% for all temperatures. Consequently we use 
iV = 5 for all computations. Inserting Eq. (131 into the 
left hand side of Eq. ([^ , one arrives at 


9p{9t\Vkip 


P^l^k 


— {dt^n.k ‘kipp}^k^n,k) 4 ” Pp.kk^k)^n+l,k : ( 1 ^) 

Note that the term proportional to the higher order cou¬ 
pling Xn+i,k vanishes due to Eq. (14). This is at the root 


of the rapid convergence of the present expansion scheme, 
see |23| . In our computation the expansion point is cho¬ 
sen such that Kk=o = Po,k=o- Hence it is an expansion at 
the physical minimum, and we are well within the con¬ 
vergence radius of such an expansion. 

Apart from the flow equation for the effective potential, 
we also need the flow equations of the Yukawa coupling, 
and the anomalous dimensions. They are deferred to the 
Appendix [^ 
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FIG. 3. Order parameter of the chiral phase transition: = (2po,fc=o)^^^ (left panel), and \dpo,k=o/dT\ (right panel) as 

fnnctions of the temperature in the different truncations of Table |I] 


C. Initial conditions 

It is left to specify the input parameters of the effective 
model. Within the FRG approach, this is done by fixing 
the initial effective action Fa in the vacuum. The UV- 
cutoff scale is chosen to be A = 700 MeV throughout this 
work. 

Choosing a larger A requires the inclusion of dynamical 
hadronisation as well as that of gluonic fluctuations: the 
validity range of the current model can be accessed in the 
vacuum by using the results of the QCD flows in [2S] 
at a given cutoff scale k = A as initial effective action of 
the model at hand. The results of such a computation 
can then be compared with the results of the QCD flows. 
The larger the initial cutoff scale A is, the larger the 
deviations grow. This limits the reliability of the model 
predictions. 

In turn, choosing a smaller A limits our thermal range, 
as the physics of temperatures with T/A > 1/7 — 1/5 
is severely contaminated by cutoff effects, see [2l]. This 
entails that the maximal temperature accessible within 
the current model is Tmax = 100 — 140 MeV. In the 
present work we also amend the model with temperature- 
dependent initial conditions which gives access to higher 
temperatures. This is detailed in the next Section [III D[ 

At the initial scale the effective potential is well ap¬ 
proximated by a classical potential, to wit 


with the relevant coupling Aa = A 2 ,a- For the sake of 
simplicity we have chosen the mass a = 0. This 
fixes the a-mass and is respsonsible for its minor vari¬ 
ation for the different truncations as well as a relatively 
large meson scattering coupling Aa at the initial scale, 
see Table |TT| Heuristically, if aiming at maximising the 
effective thermal and chemical potential range one has 
to utilise the meson mass parameter in order to minimise 
Aa. Such a choice leads to small mesonic contributions 
at large cutoff scales /c —>• A and hence at large T and fi. 
This simulates the melting of the mesons at large scales. 
Still, it rather is a feature of the model instead of built-in 
physics. In the present framework the melting of mesons 
is naturally induced by dynamical hadronisation and will 
be considered elsewhere. Here, we stick to (16). 

The other two relevant couplings are the Yukawa cou¬ 
pling h\ and the coefficient of the linear breaking term 
ca. The three relevant couplings of the model are de¬ 
termined by fitting hadronic observables in the vaccum, 
the TT decay constant /^ = ct with /^r = 93.5 MeV, 
the TT-meson mass = 135 MeV, and the quark mass 
rUg = ll2ha with mg = 303.5 MeV. The three couplings 
at the initial scale k = A as well as the predicted cr-meson 
mass at vanishing cutoff fc = 0 are summarised in Table |TT| 
for the different truncations discussed in Section imi 


D. Pressure and entropy 


VA{p) = fp\ (16) 


Truncations 

Aa 

flA 

> 

O 

1 

O 

X 

rria (MeV) 

LPA 

43 

6.5 

1.7 

585.3 

LPA' 

72 

6.5 

1.98 

584.7 

LPA + hk 

36.8 

5.5 

1.7 

551.2 

full 

77 

7.1 

2.0 

578.9 


TABLE II. Input parameters and the predicted cr-meson mass 
with different truncations. 


The thermodynamical potential density, i.e. grand po¬ 
tential density, Hfc(p) is related to the effective potential 
through 

^kip) = ^k,T[p] — r/c_T=o[p] , (17) 

which is normalised to zero at vanishing temperature. It 
has been argued above that the thermal range is given by 
A/Tmax « 5 —7, |24j. The ultraviolet cutoff A used in the 
present work has been chosen as A = 700 MeV in order 
to minimise the systematic errors of missing glue fluctu¬ 
ations and additional four-fermi and other higher order 
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FIG. 4. Pressure (left panel) and entropy density (right panel), normalised to their Stefan-Boltzmann values, as functions of 

T. 


couplings that can be described within dynamical hadro- 
nisation. However, then temperature-dependent initial 
conditions are required for T > 100 — 140 MeV: Already 
the flow at the initial scale is sensitive to thermal fluctu¬ 
ations for these temperatures, and hence the initial ac¬ 
tion as well as thermodynamical quantities, such as the 
pressure, entropy, fluctuations etc. are. This has been 
discussed at length in [22], where it has been used for a 
systematic error estimate. In the present case the ther¬ 
modynamical potential density at A can be estimated by 
the integrated flow of the difference in ( 1^ without feed¬ 
ing back the flow of the parameters. This leads to 


Ha 


4NcNf 

127r2 



1 

.expf^)+T 



(18) 


where we have dropped the negligible contribution of the 
mesonic degrees of freedom for cutoff scales above the 
chiral symmetry breaking scale <C A. We estimate 
that this procedure effectively extends the thermal range 
to T < A/5 — A/3, as it mimics the thermal change of 
the initial condition aX k = A. For the ultraviolet cutoff 
A = 700 MeV used here this gives us access to temper¬ 
atures with « 140 — 230 MeV. The corresponding 
systematic error estimates from [21] will be used later for 
an error estimate of the bayronic fluctuations. 

The pressure and entropy density follow readily as 

dp 

p = -^k=o{Po,k=o), and s = — , (19) 

from the grand potential H. Fig. (left panel) com¬ 
pares the temperature dependence of the order parame¬ 
ter /^ = (t on T in different truncations. One observes 
that the truncations cluster into two groups: full and 
LPA' in one group, LPA and LPA + h^ in the other. The 
crossovers in LPA and LPA -|- hk are broader than that 


of full approximation and LPA/ This entails that the 
non-trivial momentum dependence of the propagators, 
encoded in the wave function renormalization factors, 
increases the strength of the crossover and lowers the 
pseudo-critical temperature. This is seen very clearly in 
the right panel of Fig. where \dpQ^k=o/dT\ as a function 
of temperature is depicted. We also expect larger changes 
in LPA^ and full approximation for the high-order fluc¬ 
tuations, see Section [III E[ 

In Fig. j^we show the pressure and entropy density as 
functions of the temperature in different truncations. We 
have normalised them with those of massless ideal gas of 
quarks, i.e. the Stefan-Boltzmann (SB) limit. 


= 2iVciVy 


Ttt^ 

1 

/■ p-\ 

" 1 ^ 

(py] 

360 

12 


^ 247r2 

[tJ 


■ ( 20 ) 


One observes that the pressures calculated with trunca¬ 
tions full and LPA^ increase with T more rapidly than 
those with the other two truncations, and their values 
are also larger. These findings indicate that propagators 
with improved momentum-dependence stiffen the equa¬ 
tion of state. A similar behavior is found in the calcula¬ 
tions of entropy density, as shown in the right panel of 
Fig. ID The entropy obtained with wave function renor¬ 
malization factors rises more rapidly. 

In approximations with scale-dependent wave function 
renormalisations both the entropy as well as quadratic 
baryonic fluctuation show a peak at about 170 -180 MeV. 
This is an artefact of the current approximation, and can 
be traced back to the systematic errors discussed before. 
First of all these temperatures are above the estimate 
from [21] for the maximal temperature « 100 — 140 
MeV as well as in the border regime of the enhanced esti¬ 
mate with « 140 — 230 MeV. Hence they are subject 
to cutoff effects. Moreover, we already know from QCD 
flows at vanishing temperature that the present local ap¬ 
proximation with scale-dependent wave function renor¬ 
malisations without full momentum dependence overes¬ 
timates the anomalous dimensions in the scale regime 
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T [MeV] 

FIG. 5. Quadratic (left panel) and quartic (right panel) baryon 
different truncations, see Table |T] 



number fluctuations as functions of T in the QM model for the 



T [MeV] 


FIG. 6. Kurtosis of the baryon number distribution as a func¬ 
tion of the temperature in the QM model. 


slightly above the chiral scale, see [5^ and [521 ^ 

comparison of the scale and momentum-dependence. A 
related comparison concerning the scale dependent run¬ 
ning strong coupling with and without momentum de¬ 
pendence can also be found in [35], see Fig. 12. This calls 
for a refined analysis with momentum and frequency de¬ 
pendent propagators which will be considered elsewhere. 


Finally, we remark that the LPA^ and full approxi¬ 
mation leads to a very small negative pressure for small 


temperatures T < 80 MeV, see the inlay plot in Fig. 15 


Its origin can be traced back to the anomalous dimen¬ 
sion of the qaurk, and hints at a -minor- overestimation 
of thermal flucutation effects of the quark propagation 
for constant anomalous dimensions. We expect this ef¬ 
fect to disappear once the full frequency and momentum 
dependence of the propagator is taken into account. In 
any case, it is a very minor effect, as shown in Fig. |I5[ 
and only exists the low-temperature regime far below the 
pseudo-critical temperature. 


E. Baryon number fluctuations 

The baryon number fluctuations are usually described 
by the generalized susceptibilities, 

yB ^ A (21) 

with baryon number chemical potential /tb = 3/i. Xn’s 
are related with the cumulants of baryon multiplicity dis¬ 
tributions. Denoting the ensemble average of the baryon 
number with (Nb), the quadratic and quartic fluctua¬ 
tions are given by 

X? = yp ((“b‘) - 3(ilVB“)") , (22) 

where V is the volume of the system, and 6Nb := Nb — 
(Nb). In fact, {SNb'^) is just the width of the Gaussian 
distribution, denoted specifically by The kurtosis of 
non-Gaussian distributions k (the skewness is vanishing 
at liB = 0) is defined by 


Ka 


2 


X2 


(23) 


The quadratic baryon number fluctuation yf as a func¬ 
tion of T with different truncations is presented in left 
panel of Fig. As expected, yf in LPA and LPA -|- hk 
signal a broad crossover, while in the other two trunca¬ 
tions show a steeper transition. In general, the yB behave 
very similar to s/ssb shown in Fig.|^ 

So far, we have compared the pressure and its low- 
order derivatives with respect to the temperature (s: first 
order) and chemical potential (yf: second order) within 
different truncations. All the comparisons indicate that 
the nontrivial momentum dependence of propagators en¬ 
hances the transition strength of the chiral crossover. 

















This effect should be even more pronounced in higher- 
order thermodynamical quantities, such as x®. This ex¬ 
pectation is confirmed in Fig. right panel and Fig. 
where xf the kurtosis /tcr^are shown. Both y® and 
Ka^, calculated in the full and LPA^ approximations, vary 
rapidly during the crossover. In contradistinction, those 
obtained from LPA and LPA -|- hk as functions of T are 
much smoother. Note that the difference is quite large in 
the quartic fluctuations, in contrast to that observed in 
low-order thermodynamical quantities. 

In light of an experimental observation of (critical) 
fluctuations this is good news. It has been shown pre¬ 
viously, that effects of mesonic multi-scattering in the 
medium, as included in the FRG approach in a LPA trun¬ 
cation, soften the chiral crossover and suppress thermal 
fluctuations compared with the mean-field approxima¬ 
tion, see e.g. [28] for the kurtosis. This limits the chances 
of an experimental observation of critical fluctuations and 
hence the critical point. However, our present calcula¬ 
tions reveal that the nontrivial quantum and thermal dis¬ 
persions of quarks and mesons have the opposite effect: 
they increase the transition strength of the crossover and 
enhances thermal fluctuations. While this effect was al¬ 
ready visible in [53 [51] for the order parameter, it is 
particularly obvious for high-order thermal fluctuations, 
such as non-Gaussian fluctuations of the baryon number, 
see Fig. [^ 


IV. QCD-ENHANCED 
POLYAKOV-QUARK-MESON MODEL 

The quark-meson model studied so far is a good and 
clean testing ground for studying the relative strength 
of flcutuation effects such as multi-meson scattering and 
non-trivial dispersions. However, it does not include 
the confinement dynamics that is necessary to obtain 
the correct physics of baryon number fluctuations in the 
hadronic regime at temperatures below Tc- 

As already argued in Section|ll]the quark-meson model 
can be embedded in low energy/low temperature QCD 
by adding a background gluonic potential, see (|^. This 
directly leads to the usual Polyakov-loop-extended chi¬ 
ral models, such as the Polyakov-Nambu-Jona-Lasinio 
model [5DH55] and the Polyakov-quark-meson (PQM) 
model [53] • As explained in Section |h[ at sufficiently low 
scales the gluon fluctuations decouple due to the QCD 
mass gap, and the gluon field can be treated as a tempo¬ 
ral background field. Its contribution to the thermody¬ 
namics, however, is prominent, via the expectation value 
of the traced Polyakov loop, viz. 

Lix) = ^ {TTVix)) , and Z = ^ (Tr , 

(24) 

with 

V{x) = V exp (j,g J drAo(a;, r)^. (25) 


The expectation value of the traced Polyakov loop, L, L 
serve as order parameters for the deconfinement phase 
transition. 

In the presence of a gluonic background field the quark 
contribution to the flow is changed. This is accounted for 
by simply replacing the fermion distribution functions 
np^x^T) = l/(exp(a:/T) + 1) in Appendix [A] with the 
Polyakov-loop modified ones, i.e. 


np{x, T, L, L) = 


l + 2Le^/^ + 

1 -l- 3L -l- 3L + g3a;/T 


(26) 


In (261 we have introduced the notation 


Zq 


a: = — (I+ (27) 

Zq 

with Zq = 1, see Appendix [A] The distribution function 
of the anti-quark, 71 ^( 2 :, T) Ganges to 

np{x, T) —>■ np(x, T, L, L). (28) 

In this work we present results in two different approx¬ 
imations. In the first approximation we aim at exploit¬ 
ing the existing results of the Polyakov loop expectation 
value, e.g. from the lattice or continuum QCD computa¬ 
tions. Then the Polyakov loop expectation value serve as 
a background in which the flow of the matter fluctuations 
is computed. Furthermore, L = L is assumed and their 
dependence on the small chemical potential is neglected. 
Here we simply test the consistency of such an approach 
as well as evaluating the sensitive input parameters such 
as the field derivatives of the Polyakov loop potential. 
The respective results in this approach are collected in 
Appendix [^ 

The second approximation utilises a given Polyakov 
loop potential, and the equations of motion for L and 
L are solved together with that of the cr-meson. The 
matter fluctuations are computed on this self-consistent 
background, leading to a fully coupled system. This ap¬ 
proach is numerically more challenging, in particular if it 
comes to extensions of the present approacimation. How¬ 
ever, it has the advantage of self-consistency which turns 
out to be of eminent importance for the computation of 
baryonic fluctuations, or higher moments in general. 


A. Glue potential and higher moments 

Now we introduce the simple polynomial Polyakov loop 
model potential U{L,L), see [5T]. We shall also discuss 
the model parameters relevant for a reliable extraction of 
higher moments in QCD. A more thorough investigation 
of the model dependence of the results will be consid¬ 
ered elsewhere. The dimensionless Polyakov-loop poten¬ 
tial Vym [L,L;t] is obtained by a rescaling oiU{L,L) with 




It can be expressed in terms of the dimensionless 
reduced temperature t = {T — Tcr)/Tcr and reads 


9 


Vym{L, L; t) — 


U{L,L) 

rpH 

2 6 ^ 




(29) 


Only the coefficient & 2 (i) of the quadratic term depends 
on temperature. For large temperatures we have &2 > 0, 
while for small ones we have 62 < 0. Its temperature de¬ 
pendence above is fixed by the Yang-Mills expectation 
value of the Polyakov loop computed on the lattice and 
the Yang-Mills pressure. It has a divergence for t —>■ —1, 
that is for vanishing temperature T —0. 


62 (t) = ao 


Oi 


02 


03 


l + t (l-ft )2 (l-kt)3 


(30) 


with the parameters oq = 6.75, oi = —1.95, 02 = 2.625, 
03 = —7.44, 63 = 0.75, 64 = 7.5. Our computations 
here are based on a QCD-enhanced glue potential, that 
has the correct thermal dependence on the reduces tem¬ 
perature of QCD instead that of Yang-Mills theory. It 
has been shown in naimiii] that in the vicinity of the 
phase transition the thermal QCD scaling is obtained 
from that in Yang-Mills theory in (291 through a simple 
linear rescaling of t, 


^YM(^glue) ~ 0.57 tglue ■ (^1) 

Thus, the QCD-enhanced glue potential reads 

Fglue(T, T; tgiue) — bvM (L,L;0.57tgiue), (32) 

with tgiue = (T-TS'“)/TS'“ and « 200-210MeV 
for two flavours [22] . Here we choose 208 MeV as already 
used in [55] . 

The QCD enhancement of the potential provides the 
correct temperature scaling of the glue potential in QCD. 
QCD-enhanced computations of thermodynamical quan¬ 
tities, i.e. pressure and trace anomaly, agree remark¬ 
ably well with recent results from lattice QCD for the 
Nf = 2-1-1 flavour case already in LPA, see [22]. More¬ 
over, despite the absence of the strange quark fluctu¬ 
ations, already the Nj = 2 flavour thermodynamical 
observables are close to that in the Nf = 2-1-1 case, 
after the overall temperature scale is normalised, for a 
comparison of Nf = 2 + 1 to Nf = 2-I-1-I-1 see |M] . 
This supports the interpretation that the strange quark 
(and even more so the charm quark for that matter) con¬ 
tribute to the UV fluctuations above their mass scale and 
hence predominantly only change the overall dynamical 
scale Aqcd- In turn, at low momentum and temperature 
scales the strange quark decouples and its fluctuations 
only leads to quantitative modifications of the thermo¬ 
dynamical observables. In the present work we shall test 
this interpretation also for the baryonic fluctuations. 



FIG. 7. Polyakov loop as a function of T. The result is ob¬ 
tained by employing the QCD-enhanced glue potential in [22j . 


The respective Polyakov loop expectation values 
within the present self-consistent two flavour computa¬ 
tion are depicted in Fig. as a function of T. We em¬ 
phasise again that the potential itself is a model poten¬ 
tial that reproduces the pure glue pressure and Polyakov 
loop expectation values. It may fall short of reproduc¬ 
ing higher Polyakov loop correlation functions. The lat¬ 
ter information can e.g. be introduced to the models by 
computing higher order correlation functions or Polyakov 
loop fluctuations, see liSIEB]. A full resolution of this 
task is given by either computing the full glue potential 
from functional methods, or by computing 

the full potential on the lattice [5011641 . The full QCD- 
embedding with the computation of the glue potential 
with the FRG will be considered elsewhere. 

In any case it is of chief interest to evaluate the sen¬ 
sitivity of the results on the higher moments of the 
baryon multiplicity distribution on the higher order L, L- 
derivatives of the glue potential. We will see in the next 
subsection that the absolute values of the higher mo¬ 
ments X®, in asymptotic T/Tc-regime, crucially depend 
on parameters of the glue-potential Ogi„e- Moreover, the 
transition regime depends on the relative strength of the 
couplings in Qgi„e and This emphasises the neces¬ 

sity to compute and utilise quantitative versions of the 
glue potential, as well as a detailed study of the input- 
dependence. This is beyond the scope of the present work 
and shall be discussed elsewhere, [55] . 

Here we put forward an analytical method for the com¬ 
putation of higher moments that clearly demonstrate the 
relation between Xn and the L, L expansion coefficients 
(higher moments) of the glue potential. The higher mo¬ 
ments Xn proportional to total /r-derivatives of D. 
Hence, it requires the computation of and d^L from 
the correlation functions, that are the field derivatives of 
n at vanishing /r. This rather technical computation is 
deferred to Appendix]^ here we only briefly describe the 
main results. 
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Total /^-derivatives can be written as partial fi- 
derivatives at fixed fields and field-derivatives multiplied 


with with $ = {L,L,q,q,a,Tr), see (|D9|) in Ap- 


pendix[D| 




(33) 


Here Gij stands for the propagator, see (D7) and 


stands for a mixed partial /r and $i-derivative, see (D7|. 


Ul 


The indices i,j sum over species of fields, momentum or 
space-time indices as well as internal indices. Within this 
condensed notation, the second moment takes the form 


(D 8 ) 






(34) 


In the low temperature limit both terms, r 2„2 as well as 


that proportional to contribute. We re-order (34) 
and concentrate on xf. Note also that in the present 
approximation to the full theory we have f 2 ^,r = = 0 . 

This leads us to 




— — . 
9 


(35) 


with 


AH 2 —2 {yifia^fiLGcrL + 


+ n^^^G^L + ^iiG 


LL + 


LL ■ 


(36) 


The term An 2 carries the direct sensitivity to the model 
input, the choice of the Polyakov loop potential = 
T'^Pgiue- Most importantly this is reflected in the prop¬ 
agators Gll,Gii,G]^i that depend in leading order on 
f^giue- Note that also the other propagators depend on 
rigiue as nij is not diagonal, but these dependencies 
are sub-leading. In summary, (36) entails that for an 


accurate determination of the second moment X 2 
glue potential used in model computations should re¬ 
produce the (connected) two-point correlation functions 
Gll = {LL)c, Gj^i^ = {LL)c, and Gq^i = {LL)^ in QCD. 
Alternatively the subdominance of these contributions 
has to be shown. 

The equation for X 4 or rather in terms of 

partial /t- and 4>i-derivatives is given in Appendix in 
(D13|. Similarily to (35) it can be split in all terms with¬ 


out L-derivatives of 17 and the rest, AI74. Apart from the 
dependence on the two-point correlation functions Ar 24 
also carries a dependence of the three- and four point 
correlation function of the glue potential. 


B. Confinement and the hadronic phase 


Polyakov loop potential necessitates a detailed discus¬ 
sions of the mechanisms and scales behind the asymptotic 
behaviour of higher moments in the asymptotic regimes 
T/Tc —>■ 0,oo. A particularly interesting regime, how¬ 
ever, is the crossover regime at about as the fluctu¬ 
ations there may give indirect access to the question of 
the existence of a critical endpoint in the phase diagram. 
It is therefore very important to understand, the quanti¬ 
tative knowledge of which couplings is required to make 
predictions there. It is also of chief interest to pin down 
the relevant mechanisms in view of the model input, that 
is not fully controlled. This gives us both, access to the 
systematic error of the computation as well as the physics 
in the transition regime. 

A particularly helpful observables within this discus¬ 
sion is the kurtosis xf/xf- T/Tc —)■ 0 it simply 

counts the hadronic degrees of freedom. Hence we expect 
X 4 /X 2 1- III turn, for T/Tc —)■ 00 the kurtosis tends 

towards its perturbative behaviour, X 4 /X 2 2 /37r^. 

While the high temperature dependence is easily un¬ 
derstood in terms of perturbation theory, for T/T^ —)■ 0 
the common picture is the following: for these tempera¬ 
tures the Polyakov loop expectation values tend towards 
zero, L, L —> 0. Accordingly, the quark thermal distribu¬ 
tion np{x, T, L, L) in the presence of the confining back- 
round signaled by L, L « 0 supposedly tends towards a 
baryonic one. Indeed, taking L, L = 0 in (26) leads us to 




1 -|- 3L -|- 3L -|- 


I 


oSxjT 


(37) 


The right hand side of (37) simply is a baryonic distribu¬ 


tion function. Hence, a vanishing Polyakov loop expecta¬ 
tion value would effectively lead to baryonic properties. 
At /i = 0 we deduce from (27), 


X = —a/I 


%.fc ^ rUg 


with Too con = X - 


ha 

2z, 


fc =0 

(38) 


is bounded from below by the sum of the temperature- 
dependent constituent quark mass fc = 0. The 

—0 limit in (37) then entails, that —>■ 00 and 

-SxjT quark part of the grand potential 17 

npix, T, 0, 0) 


T 


np —t e 

is proportional to a sum of np{x, T, 0,0) 
and the 2nth-order /(-derivatives simply pull down 3^", 
leading to the desired result of counting baryonic degrees 
of freedom. Indeed this is seen when directly using L = 

Appendix 


19 


L = 0 in the model, see Fig. 

Note however, that the above picture involves taking 
L,L ^ 0 ior T —0 and then taking to 00 . A close 
inspection reveals that this explanation falls short: both 
l imi ts relate to T —>■ 0 and the thermal distribution in 
(37) depends on products of L,L and powers of 


Indeed, a necessary ingredient for the validity of (37) is 


Le^ 


JT 


< 1 . 


(39) 


The analysis of the last section concerning the sensi¬ 
tivity of the fluctuations observables on the details of the 


We shall argue now, that (39) does not hold anywhere, 
see Fig. [^for the PQM model under investigation here. 
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T [MeV] 


FIG. 8. L as a function of temperature, obtained 

in the full computation. 


We emphasise that the following argument is done with¬ 
out explicit reference to the specific model and also holds 
in QCD. To that end we discuss the equations of motion, 


= 0 . (40) 


In a flow equation approach to QCD, comprises 

the contributions from the gluonic and ghost loop in 
Fig.[TJ while comprises that from the quark and 

mesonic loop in Fig. [l] The representation in terms of 
a flow equation facilitates the discussion of the respec¬ 
tive scales. Dropping the grand potential reduces 

to the pure glue potential with a transition temperature 
Fgiua ~ 200 MeV, see [HI |2T1122]. The glue potential (in 
the current approximation) features a first order phase 
transition, and below Tgi^e the Polyakov loop expectation 
vanishes. The possible minor center-breaking corrections 
of the glue propagation due to the vacuum polarisation 
are neglected. We shall see, that for the present argu¬ 
ment we can safely drop these contributions as they only 
increase the expectation value of the Polyakov loop. 

We conclude that the non-vanishing expectation values 
of L and L in QCD are dominantly triggered by the quark 
loop in Starting at T = 0, where the Polyakov 

loop vanishes, we can estimate its T-dependence by com¬ 
puting the thermal derivative. To that end we use the 
dimensionless potential 

V[$BoM] = ^ = Klue+V_tt, (41) 


where we have extracted the canonical scaling. With 
(411 the T-derivative of (40) normalised with 1/T^ reads 
for small T/T^ —?> 0 


we are led to 


BtL = — 


-I- Vtl 
Vll + Vll 


(43) 


Now we explore the existing knowledge in QCD about 
the curvature coefficients that contribute to (43) as well 
as the mixed derivative Vtl'- 

First of all, in the limit T/Tc —>■ 0 the glue poten¬ 
tial 14 i„o,r settles to a finite potential with non¬ 

vanishing curvature coefficients Vll, V^i, Vn, see e.g. 


[58]. Hence the denominator in (43) settles at a finite 
temperature-independent value. It is left to provide a 
lower bound for the temperature-dependence of the nu¬ 
merator. The temperature-dependence of the expecta¬ 
tion values of a decays at least with the smallest mass 
scale in QCD, the pion mass. For example, in the PQM 
model considered in the present work, the only non¬ 
vanishing expectation value apart from L, L is that of 
the sigma meson, ctq. In turn, the mixed field deriva¬ 
tives, e.g. VLa can only couple via the quark sector and 
hence decay exponentially at least with 

It is left to give an estimate for the last term in the 
numerator, Vtl- To that end we split the potential V in 
a center- symmetric part and a center-breaking one. The 
lowest term in L, L in the center-symmetric part is pro¬ 
portional to LL and its L-derivative is proportional to L. 
In turn, the center-breaking terms originate in the quark 
dynamics and their contributions decay exponentially at 
least with as is the case in the present model. 

In QCD the fermionic loop in the model is substituted 
by a fully interacting one and the constituent quark mass 
in the decay goes to half of the free energy energy of an 
interacting quark-anti-quark pair. This leads to the final 
estimate for the temperature-dependence of L at /r = 0, 


OtL > max e ^ Cl 


(44) 


with prefactors Ci and Cl that carry a potential polyno¬ 
mial dependence on temperature. As the Polyakov loop 
expectation value is vanishing at T = 0 we deduce from 
(44)that 


lim L —>■ c(T) e 

T/Ta^O 


JT 


(45) 


where the prefactor c{T) carries a polynomial tempera¬ 
ture dependence. Consequently we do not reach (39) but 
rather have 


> 1 , 


(46) 


which concludes our argument. We emphasise again that 
we have not resorted to any model properties but only 
have used the scaling properties of the glue and matter 
parts, and VA.tt = respectively, 

of the grand potential. 


BtBlV = VllBtL -I- VtlBtL -I- VLit>idT^i + Vtl = 0, 

(42) 


C. Confinement and baryonic fiuctuations 


where Vtl = and index i runs over all fields 

without L, L. At /X = 0 we have L = L, BtL = BtL and 


The analysis in the last section entails the failure of the 
standard statistical confinement picture based on (37). 
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T [MeV] 

FIG. 9. Ratios of np{mq^^aa,T,L,L)/np{mq^c,a„,T,0,0) 
(full over baryonic distribution) and 

np{mq^caii,T, L, L)/np{mq^cari,T,l,l) (full over quark 
distribution) as a function of temperature, obtained in the 
full computation. 


Indeed, the latter limit does not hold for any temper¬ 
ature and the Polyakov loop augmented thermal quark 
distribution never comes even close the baryonic one, see 
Fig.|9l 

Instead one can deduce, that the hadronic limit is 


reached if (46) is satisfied. In the PQM model used here 
this happens for T/Tc « 0 . 8 , below which we only ex¬ 
pect to see the trivial thermal scaling. It is left to show 
that in this limit we see the hadronic nature of the low 
temperature phase, that is in particular X 4 /X 2 1 : ‘i®- 

spite the the thermal distribution for T/Tc settling well 
in between the quark and baryon distribution, see Fig. 
Indeed, at /r = 0 and T/Tc —t 0, and with (45), the 


thermal distribution at fc = 0 is well approximated by 

nFix,T,L,L) ^ ( 47 ^ 


potential and introduces an explicit breaking of center 
symmetry. The exponential decay factors such as 
would carry free energies of half of an interacting quark- 
anti-quark pair at infinite distance. Here we simply ne¬ 
glect the binding energy and take that of a free quark- 
anti-quark pair with constituent quark masses as a lower 
bound. This is also the approximation underlying the 
current model approximation. As we shall see, the pre¬ 
cise values of the parameters fo 2 ,& 3 ,ci are irrelevant for 
the T —)■ 0 asymptotics of the kurtosis. We also have 
dropped the mesonic terms as subleading. In ( [4^ , T is 
absorbed into x and x for brevity, and L, L on the EoM 
read 


1 


T=-\cx (e-® + 2 e- 2 ") - hL^] + 0 (e 
02 


-3- 




L=-[ci (e-" + 2 e- 2 ®) - b^L'^] + 0{e 
02 


-3- 


(49) 


Eq. (49) makes explicit the expansion in powers of 


exp{—TOq,co„/T}. All terms dropped in (48) are at least 


of the order exp{—4mq_„„/T}. It also supports the nice 
heuristic interpretation of the Polyakov loop expecta¬ 
tion values L,L being related to the free energy of a 
test quark/anti-quark: the leading order of L,L decay 
exponentially with exp{—(TOq^,,„„ ± /r)/T} respectively. 
The implicit and terms encode terms of at least 
quadratic order in exp{—(wg^con ± Eq. (49) also 


entails that at finite chemical potential L is enhanced rel¬ 
ative to L. Inserting the respective equations for L, L on 
the right hand sides of (49) leads us to 


r —a 

02 


r —a 

O 2 




where x in (j^, (38) is bounded by mq^^^^/T for all k. 
Hence the tnermal distribution decays with a thermal 
correlation length related to the inverse of twice the con¬ 
stituent quark mass. Note also that the quark contri¬ 
bution of (the flow of) t he g rand potential H is directly 
proportional to np, see (A5) in Appendix [A| 

With the scaling properties for T/Tq —)■ 0 deduced 
above we can expand the grand potential in orders of 
exp{— 77 iq^eon/^r} on the equations of motion. This leaves 
us with 


Eq. (50) reveals the composite nature of the Polyakov 


loop expectation values. For example, L has an expan¬ 
sion in quarks and anti-diquarks all of which carry the 
same SU{N) representation. Hence we expect LL to 
carry both mesonic {qq) and baryonic properties {qqq, 
qqq). Indeed, it follows from (50) that 


LZ=^e-2- 

^2 


H = - —LL - — (L^ 
2 6 ^ 


+ Z(e 


2 e- 


Cl 

2 


+ flmes + O 


L (e 


2 e 


-2x\ 


(48) 


where the 62 , ^3 terms stem from the expansion of the 
glue potential Hgiue in powers of A, Z, see e.g. [53], and 
we have dropped the subleading term 64 / 4 (LZ)^. The 
ci-term stems from the quark contribution to the grand 




—3x 


+ e-^4 


(51) 


The first term in (51) carries no baryon number and 


heuristically relates to mesons, while the second one car¬ 
ries baryon number, there are no diquark contributions. 
The occurance of the decay with multiples of the con¬ 
stituent quark mass for qq and for qqq relates 
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to the fact that the Polyakov loop terms describe pairs 
and triplets of quarks rather than the bound states. 

In case of the baryons this is related to the missing 
higher order quark interactions in the current model that 
carry the interaction energy of the hadronic states. This 
can be amended by taking into account multi-quark in¬ 
teractions with or without dynamical hadronisation. 

In case of the mesons the missing contributions are 
hidden in the mesonic part of the grand potential. To 
see this more clearly we proceed with the expansion of 
the full grand potential. Note first that 11 in (48) only 


has linear and bilinear terms in L and L evaluated on 
their EoMs, ( |4^ , as well as cubic terms. The latter cu¬ 
bic terms carry baryon number in leading order. Using 
[LBl + Ldi)ft = 0 to express the linear terms with the 
bilinear and cubic ones we arrive at 


n =^LL + y (i' + L^) + + o 



T [MeV] 


FIG. 10. Pion decay constant as a function of T, calculated 
in the PQM with LPA and full approximation. 


— ^nnn “t“ Hflo Hn 


(52) baryons. Eq. (54) implies in particular 


with 


^qqq — 




—3x\ 


o 




1 ^- 2 - 


n - ——p 

liqq - 2 ^^ e 


0 {e 


-4- 


4- 


(53) 


Note that the EoM also imply that the linear terms have 
the for m in (51) with additional contributions in the first 
line in (531 steming from the +L^-term. With (52) the 


grand potential takes its final baryonic-type form. It has 
a /r-independent offset and subleading contributions. The 
flqqq term carries baryon number with contributions from 
all three terms in the grand potential. In terms of the 
Polyakov loop it has the form + L^. The mechanism 
unraveled above only holds for L L —> 0 which signals the 
hadronic phase. 

Note also that the total mesonic contribution is given 
by + The latter part dominates in the hadronic 
phase and shows the correct mesonic mass scales. We 
expect a similar phenomenon for the baryonic part once 
the baryonic states are included in the effective action. 
Note however, that the baryons do not affect the dy¬ 
namics of the theory due to their heavy masses which 
suppressed any baryonic off-shell fluctuations at all mo¬ 
mentum scales. 

Now we are in the position to discuss the fate of the 
kurtosis deep in the hadronic regime, that is for low tem¬ 
peratures T/Tc «C I. Prom (52) we deduce, that even 


powers of total /r-derivatives of the grand potential are 
simply proportional to the potential in the T —> 0 limit, 




.-4- 


(54) 


-^ = l + 0{e 
X2 


(55) 


In summary, despite the final outcome being the same, 
the mechanism behind the asymptotic T-scaling of the 
kurtosis and also the higher moments is different to that 
suggested in the literature: in contradistinction to former 
claims the asymptotic low energy scaling is obtained in a 
regime where the thermal distribution are far from being 
baryonic, see Figs. |8|9[ Indeed, the Polyakov loop en¬ 
hanced thermal distribution is far from the baryonic one 
for any temperature, as is clear from the above analysis. 


D. Results 

In the last section we have reached a better under¬ 
standing of the mechnisms taking place in the hadronic 
low energy/low temperature regime of QCD. Most of the 
properties discussed do indeed not depend on the model 
under investigation but are properties of QCD. With this 
in mind we hnally come to the numerical computations 
of this paper for the QCD-enhanced Nf = 2 PQM model. 
Before we present the results, let us briefly summarise the 
approximations used here: We build upon the approxi- 


with flqqq counts three-quark states, the effective 


mation tested in Section III for the quark-meson model. 
We use a full effective meson potential I4((/>;, L, L) that 
includes multi-meson scatterings. The dependence on 
the Polyakov loop variables L, L is implicit as the flow 
is evaluated on given backgrounds L,Z. We also in¬ 
troduce scale-dependent wave function renormalisations 
Zq,(L, L), Zq{L, L) with the same implicit dependence on 
L,L. The Z’s account for the change of the dispersions 
for both mesons and quarks. Additionally we allow for 
scale-dependent quark-meson interaction hk{L,L). The 
quark-meson part in the L, L background is amended 
with a QCD-enhanced Polyakov loop potential I4i„e(T, L) 
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FIG. 11. Pressure (left panel) and entropy density (right panel) as functions of T/T^. The psudo-critical temperature is 
determined by the peak of |9po,fc=o/3T|, Tc = 181 MeV, 190MeV for fnll and LPA, respectively. 



FIG. 12. Left panel: trace anomaly as a function of T. Right panel: trace anomaly for N[ = 2 + 1 in LPA from (red dashed 
line), multiplied by the ratio of the full A^/ = 2 data over the LPA data (black solid line), in comparison to lattice result from 
Wuppertal-Budapest collaboration m- 


in Eq. (321. The potential is a standard model poten¬ 
tial with the correct temperature scaling of QCD, see 
[M1211 US- Our full results are compared with results 
in the standard LPA approximation as has been used in 
[22] for Nf = 2 and Nf = 2-1-1 flavour QCD. 

In Fig. we show the pion decay constant as a func¬ 
tion of temperature. As in the QM model the additional 
fluctuations from the wave function renormalisation and 
the Yukawa coupling lead to a sharper chiral crossover. 
In [22] it has been shown that in the 2-1-1 flavour case 
there is a very good agreement of both pressure and 
trace anomaly as functions of reduced temperature for 
T/Tc < 1.2 with the corresponding lattice results of the 
Wuppertal Budapest group, see Figure 5 there. Note that 
the pressure even agrees well above this scale. However, 
the trace anomaly, which is more sensitive to the correct 
counting of fluctuations, shows a widening gap. The LPA 
results overshoot the lattice results which hints a both, 
an incorrect counting of fluctuations at large scale as well 
as the impact of the ultraviolet cutoff. 


As already mentioned before, the agreement between 
the 2 - 1-1 and 2 flavour results as functions of reduced 
temperature is still remarkable for T/T^ < 1. This sup¬ 
ports the point of view that the dominant effect of the 
strange quark on the thermodynamics at temperatures 
T < Tc is to change the overall scale. For larger mo¬ 
mentum or temperature scales, however, its dynamics 
contributes inevitably. Keeping this in mind we com¬ 
pare the LPA and full pressure and entropy in Fig. |ll| as 
functions of reduced temperature. The additional fluctu¬ 
ations introduces via the wave function renormalisations 
and the Yukawa couplings have little impact on the ther¬ 
modynamics. We expect this property to hold also for 
2-1-1 flavours. Hence the already very good agreement of 
the pressure and entropy in LPA is not spoiled. 

The trace anomaly, left panel of Fig. also does not 
change signihcantly for temperatures below Tc- Above 
Tc, the full results undershoots the LPA result. Again 
we expect the 2 - 1-1 result to show the same change. 
However, in this regime we are sensitive to cutoff effects 
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FIG. 13. Quadratic (left panel) and quartic (right panel) baryon number fluctuations as functions of T in the PQM model in 
LPA and full truncation. 


and missing fluctuations in the first place. 

We conclude that the additional quantum and thermal 
fluctuations considered in the present work have very lit¬ 
tle impact on the thermodynamics below Tc, despite their 
significant change of the matter fluctuations as shown in 
the QM model. Above they have a larger impact as 
can be seen from the trace anomaly, left panel of Fig. 

In a bold attempt to estimate the results of a 2-|-1-flavour 
computation we also show the Nf = 2-1-1 data in LPA 
for the trace anomaly from , multiplied by the ratio of 
the full Nf = 2 data over the LPA data, see right panel 
of Fig. As expected, the quantitative agreement with 
the lattice data is not spoiled for T < Tc- For tempera¬ 
tures T >Tc one clearly sees the effects of the ultraviolet 
cutoff and the missing fluctuations. 

Finally we come to the baryonic fluctuations.The nu¬ 
merical results for the second and fourth moment are 
shown in Figure We emphasise that the computation 
of the /x-derivatives of the grand potential at small tem¬ 
peratures requires the careful evaluation of the genuine 
thermal and chemical potential scaling of quark anoma¬ 
lous dimension and Yukawa coupling, for more details 
see Appendix Interestingly, the difference of the full 
results to that in LPA is significant in clear contradis¬ 
tinction to the situation for the thermodynamics. This is 
expected as the higher moments by definition carry the 
details of the fluctuation physics in the model or theory 
at hand. From Fig. one observes once more that the 
quadratic and quartic baryon number fluctuations ob¬ 
tained with the full approximation, are larger and change 
more rapidly during the chiral crossover. Note that larger 
fluctuations are important in view of experimental mea¬ 
surements, whose statistical error increase significantly 
for higher fluctuations. 

In Fig. [^we present the final results on the kurtosis. 
As expected, the significant change from the LPA results 
to the result of the full approximation propagates from 
that in the second and fourth moments X 2 wf- This 
is good news in terms of experimental measurements. 


and provides a very good test for the convergence of 
the present approximation scheme towards the full re¬ 
sult. Such an analysis also requires a discussion of the 
systematic error. The grey band in Fig. provides a 
rough estimate of the systematic error of the full com¬ 
putation arising from the UV-cutoff scale A, through the 
following simple formula |22) : 


B B B 

X 4 I A X 4 _ X 4 I 

X 2 X 2 X 2 


I± 


3A/T _ 1 


(56) 


with A = 700 MeV. Note that this error estimate does 
not include the effects of the correct hadronic decoupling 
above Tc which in the present framework is most eas¬ 
ily implemented with dynamical hadronisation. More¬ 
over, at even higher scales gluonic fluctuations have to 
be added. Finally, for temperatures T < 0.8 Tc we ex¬ 
pect frequency-dependences in particular of the fermionic 
couplings to further improve the reliability of the results. 
This estimate is based on a computation of fermionic 
couplings one the first fermionic Matsubara frequency ttT 
and the improved one used in the present work, that takes 
into account the correct thermal and /x-scaling. The lat¬ 
ter computation implements effectively the effects of full 
frequency dependences, for details see the Appendix [C| 
The results start to disagree for T < 0.8 Tc and we ex¬ 
pect a compution with full frequency-dependences to de¬ 
crease the systematic error further in this regime. A more 
detailed analysis in an upgraded approximation will be 
discussed elsewhere. 

The expected asymptotic value 1 for the kurtosis is 
obtained at low temperature in both LPA and the full 
approximation. We also compare the kurtosis of baryon 
multiplicity distribution with that of the continuum- 
extrapolated lattice result from Wuppertal-Budapest col¬ 
laboration [5]. The latter is a 2-1-1 quark flavour result 
with physical quark masses. In the spirit of the previ¬ 
ous discussion about the dominant strange quark effect 
as a different total scale setting in 2-1-1 flavour QCD in 
comparison to 2 flavour QCD we rescale the temperature 
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T [MeV] T/T, 

FIG. 14. Left panel: Kurtosis of baryon number distribution in the PQM model. The grey band shows an error estimate for the 
full computation, resulting from the UV-cutoff scale. Right panel: Comparison with the continuum-extrapolated lattice results 
from Wuppertal-Budapest collaboration [^. The temperature is rescaled by their corresponding pseudo-critical temperature 
Tc 


in units of the pseudo-critical temperature Tc- For the 
definition of the absolut temperature scale we choose the 
peak position in the fourth moment xfj this is eas¬ 
ily accessible in both approaches. From the right panel 
of Fig. ^ one can easily read = 193 MeV, 210 MeV 
for full and LPA, respectively. For the 2-1-1 flavor lattice 
simulation, Tc = 155 MeV is obtained from the results of 
X? 0, which is also consistent with the calculations in 
m- With such a normalisation one finds that the full 
calculation is in very good agreement with the lattice re¬ 
sult for temperature T < Tc, while the LPA is not, see 
Fig. 14 The most significant difference is the rapid con¬ 
vergence of the kurtosis in the full computation and the 
lattice result towards 1. This entails that the dominance 
of the hadronic nature for T < Tc is approached very 
quickly, which also has an impact on the rapid chemical 
freeze-out that is observed experimentally. The present 
approach allows to study the relative importance of glue, 
quark, mesonic and baryonic flcutuations which will pro¬ 
vide further insights in this intriguing question. 


V. SUMMARY AND OUTLOOK 

In this work we have studied the phase structure, ther¬ 
modynamics, and in particular higher moments of baryon 
multiplicity distribution, in low energy effective models 
with matter and glue fluctuations. The fluctuations have 
been taken into account with the functional renormalisa¬ 
tion group. This has allowed us to investigate in detail 
the impact of the matter and glue dynamics on baryon 
number fluctuations and the kurtosis of baryon number 
distribution. 

Furthermore, we have provided new insights on the 
mechanism of quark confinement in terms of effective 
degrees of freedom: the common picture in low energy 
effective models is that effective quark confinement is 


realized through an effective baryonic thermal distribu¬ 
tions, called statistical confinement. A detailed study 
unravels that the effective thermal distribution is quite 
far from the baryonic one, and the deviation grows ex¬ 
ponentially when T/Tc approaches zero. Instead of an 
effective baryonic thermal distribution we find that the 
center symmetry of the glue potential plays a key role, 
which guarantees that the degree of freedom is baryonic 
at low temperature. 

For a detailed study of the role of fluctuations several 
different approximations to the full low energy effective 
theory of 2 flavour QCD have been employed. These 
extended approximations have been compared with the 
standard LPA truncation, that allows for a very effi¬ 
cient inclusion of meson fluctuations and has been con¬ 
sidered in the literature. We find that the non-trivial 
frequency and momentum dependence of the propaga¬ 
tors strenghens the temperature-driven chiral crossover. 
In another words, the non-trivial dispersion introduced 
by quantum and thermal fluctuations increases the tran¬ 
sition strength. This effect is observed in all observables 
considered in this work: the order parameter of phase 
transition, the pressure, the entropy, trace anomaly and 
baryonic fluctuation observables yf and xf. Interest¬ 
ingly, while the effects on order parameter and thermo¬ 
dynamical observables is very small, it is significant for 
the higher moments. This is very well compatible with 
the fact, that the results for 2-1-1 QCD for order param¬ 
eter and thermodynamics including the trace anamaly 
were already in quantitative agreement with the lattice 
results for T <Tc for the LPA approximation. 

We have also provided estimates for the trace anomaly, 
see Fig. \V2\ and the kurtosis, see Fig. in 2-|-l flavour 
QCD. This is based on the comparison of 2 flavour with 
2-1-1 flavour computations in [22] as well as 2-1-1 with 
2-I-1-I-1 flavour computations in [^. This suggests that 
the predominant effect of strange and charm quark is 
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the simple rescaling of the total momentum scale which 
is done by presenting the results as functions of the re¬ 
duced temperature. In the validity regime of the present 
truncation to low energy QCD, for temperatures T < Tc, 
the present results agree very well with the corresponding 
lattice results. This provides further support for the reli¬ 
ability and the convergence of the present approximation 
scheme of the QCD-enhanced model for these tempera¬ 
ture regime. 

Finally we discuss our findings in view of the experi¬ 
ments: the present analysis predicts larger baryonic fluc¬ 
tuations in QCD at finite density in comparison to model 
results in the literature obtained in LPA. This effect is 
quite prominent and, thus, it may play an important role 
in the observed kurtosis of the net-proton multiplicity 
distributions at RHIC. 

Technically, it has been shown in the literature that 
meson fluctuations beyond the mean-field approximation 
smoothen the chiral crossover. This also leads to a signifi¬ 
cantly reduced critical region in comparison to mean field 
results as well as small baryonic fluctuations. In sum¬ 
mary the mesonic fluctuation analysis made the search 
for a critical end point far more difficult. The present 
analysis emphasises the relevance of fermionic fluctua¬ 
tions. The latter naturally strengthen the crossover and 
in particular enhance baryonic fluctuations observables. 
This increases both, the likelihood of a critical endpoint 
as well as its observability via fluctuation measurements. 
This is in line with the observation that the measure¬ 
ments of higher moments of the net-proton distributions 
in the first phase of BES program at RHIC show that 
their energy dependence can neither be reproduced by 
transport models without critical point contributions, 
nor by a hadron resonance gas model [68]. Hence, al¬ 
though more analyses and explanations are needed to 
clarify the measurements, the QCD critical point is still 
a fascinating candidate. 

The calculations in this work have been performed in 
a QCD-enhanced low energy effective model. Its lim¬ 
itations at T > Tc as well as the missing frequency- 
dependence of the couplings have been discussed in de¬ 
tail. Ongoing extensions of the current work concern 2-1-1 
flavour QCD, the frequency-dependence relevant in par¬ 
ticular at small temperatures T < 0.8 the missing glue 
fluctuations and hadronic decoupling at T > Tc. The lat¬ 
ter is considered within dynamical hadronisation. These 
improvements entail a full embedding of the present low 
energy effective theory in continuum QCD in the FRG 
framework as done in the fQCD collaboration [B3]. We 
hope to report on related results in near future. 
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Appendix A: Regulators and threshold functions 


In the present work we use 3d flat regulators, @11171, 
for quarks and mesons. They read 


with 


Rlido, q) = Zq^kn ■ , 

rB{x) = 0(1 - a;), 

rpix) = - 1^ 0(1 - x). 


(Al) 


(A2) 


In order to simplify our expressions for the threshold 
functions, we define 




1 


Z^do + 1 + 


kTlq 


zUqo + -I- 1 -bm; 


(A3) 


q^k 


with go = do/k, p, = ^J^,|k, and go = (2^^ -b l)7rT (uq S 
Z) for fermions and 2nqTrT for bosons. Note that z = 
which are the ratios between the longitudinal 
and transverse wave function renormalisation factors, are 
chosen to be 1 throughout this paper. Eirst of all, we 
define 


^(1) z^]T,) = J . (A4) 


Summing up the Matsubara frequencies yields 

R(l){'dlq,k,Zq;T,fl) = 


2zq^/l + mlf^ 


X - npiml i., Zq; T, n) - npiml ,., Zq\T, 


(A5) 


and 

^{i)i.'fdl,k,Z4>-,T) 
1 


zT\l^+H 


= {^+nB{mlf,,z^\T)^ , (A6) 

4>,k 








18 


with the distribution functions 




exp<j + ^ I - 1 


(A7) 


and 


nF{mli.,Zq-,T,^j.) 


(A8) 


exPi ? + 


+ 1 


The threshold functions appearing in the flow equation 
for the effective potential, i.e. Eq. are given by 




and 




(l-^)^(i)«„z, = l;T,/r). (AlO) 


d- 1 

Furthermore, we also need 

T 


^{n){ml^k.Zq-T,ii) = -J2{Gq{q,mlk)T , (All) 


which are easily obtained from through 

1 d 

^(n+l) {"^q^kJ T’, /r) = o - 2 i''^q,k^ ^95 ^5 d) ’ 

^ ^'^q,k 

(A12) 

The threshold function BB( 2 , 2 ) is derived from 
BB^ 2 ,k^ fhl^,k^ z,p-,T) 

, /c , /c 


with 


^^(1,1) i^^ipa^k't ^^ly,k^ Zffj-i T) 


T 


. (A14) 

riq 

Its analytic expression, after the Matsubara frequencies 
are summed up, is given by 


1 f/I 

a72 


+ nB{fhl^,k. T)) 


(f>a-,k) 


1 


■ + 


+nB(m^^ fe,Z0;T) 




(A15) 


In the same way, by means of derivatives with respect to 
appropriate masses, threshold functions TB's in Eq. (C4) 
are easily obtained with 


rj' _^ 

FB(^i^i){ffil,.,fhl,.,Zq,z^;T,fi,po) = - - q,fhli,)Gq{q,mlj,) 


1 f t-2 

- {nBirhl,k,z<j>:T) + l) ■ 


1/2 


+ ^l,k) (ipo-^k+^i^ + ml k^G^ _ + 


(l + "il/c)^ (ipo-P-+ - (l + "i9./c)(f) 


+ nF{m k,Zq;T,-p) 


-V 1 

+ {riF{rh1,^,Zq\T,p) - l) 


(1 + K,k){^Po - d + ^ (1 + ml ,)'/")'-(! + ml,) 


(AI6) 
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The T —0 limit of (A16) is achieved by simply drop¬ 
ping the bosonic thermal distribution functions, ub 0, 
as well as reducing the fermionic ones to Heaviside theta 
functions, nF{m^ f,,Zq\T,±^) 6»(fc/zq(l-|-m^ 

The latter limit signals the silver blaze property: be¬ 
low the critical chemical potential /r = /u* no observable 
shows a dependence on the chemical potential. Note that 
for T —)■ 0 there is still a /r-dependence on pg ~ Again 
this relates to the silver blaze property as observables 
relate to pg-integrals, and the p-dependence can be re¬ 
moved by a simple shift of the integration contour below 
the first pole in the complex plane at /r = p*. Conse¬ 
quently, on the level of the threshold functiFH(i in 
(A16) there is a qualitative difference between the fi- 


dependence in the quark distribution functions np and 
that on the frequency po + i /i: only upon frequency inte¬ 
gration both ^-dependences show the characteristic ex¬ 
ponential dependence on exp{±/r/T} relevant for the sil¬ 
ver blaze property. This is investigated in detail in Ap¬ 
pendix 

Here we proceed with the relevant threshold functions 
^5(2,1 )ml ,,, Zq, z,p; T, /r,po) 


T 




' k ^ 


G^ip- 


<l^rnlk){Gq{q,m: 


'q,k 


■)y 


Zq, z^-, T, p,po), 


and 


ml ,,, Zq, Z4,; T, fj.,po) 

T 2 

= ^'^{G4p- q,ml^k)) Gq{q,fhl,,) 

riq 


(A17) 


d 

= - F, - 2 ^-^(1.1) ml ^., Zq, Z4,; T, p,po) . (A18) 

Note that the functions FB(ix,, FB( 2 ,i)^ and FB^i ^) are 
complex valued at non-zero chemical potential and, thus, 
we have to project them onto the real axis, as explained 
in the text, viz. 




(A19) 


with z,p = Zq = 1. The threshold function in the flow of 
the Yukawa coupling (C5) can be expressed as 

-(d) 


^(1,1) Vq.k, V^,k-, T, p) 


d-l 


(1 - 


Appendix B: Anomalous dimensions of the mesons 


We have discussed the flow equation for the zero-point 
correlation function, i.e. the effective potential, in Sec- 
tion |HI A| In order to investigate more interesting proper¬ 
ties embedded in the meson and quark propagators, one 
has to go beyond this order. Differentiating both sides 
of Eq. Q with respect to fields twice, one arrives at the 
flow equations for two-point functions. After perform¬ 
ing appropriate projects, one can obtain the desired flow 
equations. For example, the meson anomalous dimension 
can be obtained with 


il4>,k{Po,p) 


1 1 S^dtTk 

Z^,k ’ 

(HI) 


where a sum over the iV| — 1 pions is implied. The flow is 
evaluated at the expansion point k of the Taylor expan¬ 
sion in the mesonic field, as the right hand sides of the 
flows are evaluated at this point. Note that another pro¬ 
jection procedure for rj relates to finite differences, e.g. 
dtT^^\po = 0,p^ = k'^) — i9tr[,^^(0), as has been used in 
[701 El]- This gives better access to the global change 
of the wave function in the relevant momentum regime 
< k^. However, it has been shown in [Mj that a 
Taylor expansion about <7 = 0 works quantitatively (for a 
mesonic O(A^)-model) for momenta q^ k'^■ This origin- 

taes in the infrared regularisation of these momenta. As 
il 4 >,k is only inserted in the right hand side of flows with 
loop momenta (f < k^ we use = 0 as the evaluation 
momentum. We also choose vanishing frequency <70 = 0 
in order to keep an expansion point that has the under¬ 
lying Euclidean 0(4)-symmetry. Note, however, that the 
regulator used breaks the 0(4)-symmetry and we will re¬ 
evaluate this choice in the case of the fermionic anoma¬ 
lous dimension. 

In summary this leads us to 


V(p,k — 


= -^{^'^kiK\kk)fBB(^2,2){ml^k,ml^FT) 

+ Njil [(277,,fc - 3) J'(2) {fhl k, T, p) 

-4(77,,fc-2)J'(3)(TO2fe;r,7i)]|, (B2) 

where the threshold functions BB( 2 , 2 ) and F(n) are de¬ 
fined in [23] , and are presented in Appendix as well. 


Appendix C: Flows of the fermionic couplings 


+ (1 ~ 

with the fixed frequency po,ex 
pendix [^ 


ml. 


,k\T, p) 


(A20) 


to be specified in Ap- 


In the spirit of the derivative expansion used for the 
present class of models the flow of the fermionic couplings 
is obtained by evaluating the corresponding derivatives 
w.r.t. the fields at low frequencies and momenta. These 
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flows are given by 


Vq,k{P0,l^ 


1 1 

Zq,k{Po,P) iNcNf 



S^dtTk 

6q{-p)Sq{p) 


(Cl) 


for the anomalous dimension and 

1 


dthkipo,:^ =- 


a ANcNf 


xRe 


Tr - 


S^dtTk 

Sq{-p)Sqip)J 




(C2) 


for the Yukawa coupling. The results in [24] for the sim¬ 
ple spatial momentum dependence for ^ of the 

mesons carries over to the quarks as it originates in the 
infrared regularisation of these momenta. Hence we use 
<7 = 0. 


It is left to specify the frequency at which (Cl I and 


(C2) a re eva luated. We have already discussed briefly 
below (A16) in Appendix [A| that the Tt-dependence in 
the argument po + i p only reflects the correct physical 
behavior if keeping the full frequency dependence. Hence, 
within an evaluation of the fermionic flows at a fixed fre¬ 
quency we face subtleties that relate to the silver blaze 
problem: all couplings with fermionic legs have an ex¬ 
plicit T and /r-dependence as well as one in the frequency 
argument. At finite temperature and chemical potential 
the latter is a sum of the fermionic matsubara frequency 
and the chemical potential, po + i p = 2TrT{n + 1/2) + ip 
with n G Z. The former explicit dependence shows the 
standard thermal suppression as well as that of the pa¬ 
rameters such as the masses and couplings. This is clearly 
seen in ( |A16 ) where the explicit p and T-dependence is 
that of the fermionic distribution functions np as well as 
that of z^/q. 

We conclude that inserting po + i ^t-dependent anoma¬ 
lous dimensions in the loops leads to standard thermal 
sums, and finally to additional thermal factors np and 
1 + np related to thermal distribution functions as in 


(A16). Related results will be considered elsewhere. 

In turn, frequency-independent wave function renor¬ 
malisations are obtained by evaluating rjq^kipo + ip^p) 
at some fixed small fermionic Matsubara frequency po = 
2TrT{n+ 1/2). This introduces an artificial temperature- 
dependence. Note also that this subtlety is not related to 
the specific regulator used but rather applies to all reg¬ 
ulators. For small cutoff scales k < T this is a sublead¬ 
ing effects, for k ^ T, however, thermal effects should 
be suppressed exponentially with exp {—c{Rk)k/T} with 
a cutoff dependent coefficient c{Rk), see |72|. For the 
present 3d flat cutoff we have c{Rk) = 1/2 leading to the 
standard thermal distribution functions already for the 
anomalous dimensions. 



FIG. 15. Pressure in unit of fm“^ as a function of T with 
truncations full and LPA. The inset plot zooms in to the re¬ 
gion T = 0 ~ 100 MeV. 


The articficial temperature-dependence is sub-leading 
as it is polynomially suppressed with powers olT/k but it 
may play a qualitative role in temperature regimes where 
the observables at hand is small. Moreover, the larger the 
canonical momentum dimension is, the bigger is the ef¬ 


fect. In Fig. 15 we show the pressure of the QM model 
in the low temperature regime. There, it is exponentially 
suppressed with exp{—TO tt/T} and hence the above cut¬ 
off effects play a role. Here, pq^k and hk are ob tain ed by 
evaluating the corresponding flows, (Cl) and (C2) sim¬ 


ply on the first Matsubara frequency ttT. The discussion 
above explains the negative pressure seen in Fig. |15| for 
the QM-model in the full approximation. 

Evidently, this problem is resolved by improving the 
approximation to frequency-dependent quark anomalous 
dimensions Pq{p) and Yukawa couplings h(j)). Here we 
resort to an effective resolution that does not force us 
to crank up the approximation: on the right hand side 
of the flows both pq and h are evaluated for momenta 
(f < k^. The integrands decay proportional to if for 
small momenta f > 0 and are maximal for momenta 
f « 1. This anyway suggests an evaluation of the 
flows for pq and h f Ki k"^. The present approxima¬ 
tion assumes 0(4)-symmetry of the effective action also 
at finite fc, and hence the couplings are assumed to be 
func tions of -I- f. Thus, we will evaluate the flows 
(Cl) and (C2) at Po = ^ for A: > ttT. This resolves the 


problem of the artificial temperature dependence. 

For k < ttT the minimal momentum squared is 
bounded by the lowest Matsubara frequency and hence 
P^ > (ttT)^. This suggests an evaluation at 


pLx = + fT)HT{k/T) , 


(C3) 


with some smoothened version Op of the heavyside func¬ 
tion with 6T{k/T ^ 1) —>■ 1 and 9{k/T ^ 1) —>■ 0 reflect¬ 
ing the known thermal decay, e.g. 9t{x) = exp{—2a:}. 
Such a 9p{2k/T) overestimates the effect of the Matsub¬ 
ara sum and we take the conservative choice exp{—2/5a;} 
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FIG. 16. Quadratic (left panel) and quartic (right panel) baryon number fluctuations in the full approximation with external 
frequency po + i fi = Po,ex, denoted by full, and that with external frequencies po + i jJ, ~ i^T + i p, denoted by full. 



T [MeV] 


FIG. 17. Same as Fig. |16[ but for the kurtosis of baryon 
number distribution. 


which guarantees the decoupling for large cutoff scales 
k/T 3> 1 and does not influence the thermal behaviour 
for k/T < 1. 


The above discussion of the temperature-dependence 
applies as well to the dependence on the chemical poten¬ 
tial. An evaluation of pq at a fixed Matsubara-frequency 
destroys the relation between p and po"derivatives and in¬ 
duces an artificial /r-dependencies. Again this may play 
a role for small temperatures where /i-derivaties tend to¬ 
wards zero. Hence, in the present work we simply drop 
the su b-leading /i-dependencies and set + i p ^ in 
(A16) and the following equations. In summary we arrive 


at the anomalous dimension 

X - l)J'6(i,2)(TOq ^,771^ ^,;T,/i,po,ex) 

+ .Po.ex)| • (C4) 

and the flow of the Yukawa coupling 

dthk =(^ + V<i,k)hk + 

x{-{N] - fc,775,fc,?70,fc;r,/r,po.ex) 

4” Vq,kj i T, J>0,ex) ^ ■ 

(C5) 


Let us now compare the results for y® j xf ^md the kur¬ 
tosis X 2 /X2 fo'' the full aproximation with consistent T 
and ^-dependence leading to (C4|, (C5) with the com¬ 
putation with fixed frequency po + i p = ttT + i p, de¬ 
noted by full. The results for the second and fourth mo¬ 
ment xf, X4 and for the kurtosis x?) X? are presented in 
Figs. [umi Evidently, both approximations are consis¬ 
tent with each other except for low temperatures, as can 
be seen in particular for the second moment in Fig. |16| 
These deviations are even more obvious in the kurtosis. 
Fig. [HI where the deviation below T « 150 MeV is sig¬ 
nificant, that is for T < 0.8 Tg. The full apropximation, 
expanded at po + * M po.ex has a T, /i-consistent scal¬ 
ing and exhibits the correct asymptotic behavior of the 
kurtosis for low temperature in Eq. (55). In turn, an ex¬ 
pansion at Po + * M po,ex + i p fails for temperatures 
T<0.8Tc. 
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Appendix D: Baryon number fluctuations and 
vertex expansions 


The generalised suceptibilities in (211 involve total nth 
/is-derivatives of the pressure p = —where $eom 
are the solutions of the equations of motion of all helds 
in the model, = (L, L, q, g, cr, tt). These derivatives hit 
the explicit /rs-dependence of the grand potential, and 
its implicit ones in the couplings as well as that of ^eom- 
In the following we discuss p = /tb/ 3-derivatives of 
the grand potential n. For example, its hrst /i-derivative 
reads 




dp 


go 

= d^\^n + (Dl) 


where the evaluatio n of the right hand side at $ = $eom 
is understood. In (Dl) we have used the equations of 


motion, and the partial derivatives are at fixed $. In 
(Dl) and the following relation we shall use, that the 


total /x-derivative can be split in a partial /x-derivative at 
fixed field and a part that hits the implicit ^-dependence 


of $eom, to wit 


dp 




^Z,EoM 


d 


With (D2) the second /x-derivative follows as 


dp^ 


= din + 


92 D 
‘a$,9/x 


(D2) 


(D3) 


Eq. (D3) depends on the 9^ <l>qEoM- This can be computed 
from the ^-derivatives of the respective EoM, to wit 


d dn 

dp d^i 




with 


(2) ^ d'^n 


This can be resolved for leading to 


^/X^2,EoM - , 


(D4) 


(D5) 


(D6) 


where, for the sake of brevity, we have introduced the 
notation 


n 


gn+m^ 




dp^dcj),^ ■ ■ ■ d(t)i^ 


and Gij = 


1 


D(2) 


(D7) 


which also keeps the higher /x-derivatives simple. Insert¬ 
ing (D6) in (D3) and using the notation (D7) finally leads 
to 


d^D[4>EoM] 

dp"^ 




(D8) 


Now we proceed to the third /x-derivative of D. Within 
our short hand notation (D7), and using (D6), the total 


/x-derivative, (D2), takes the form 

JL - f) _ o r JL 

dp~ ^ 


(D9) 


Applying ( |D9| ) to (D8) leads us to 

(DIO) 

The total //-derivative in the second term on the right 
hand side can also be performed with 


d// 




djji 


Gij — Gi]^ ^klrGrs^fis) Gij , (Dll) 


following from (D9). Then we finally arrive at 


d//3 


— 3^^jGji ^^jGji 


^fiiGir ^filGlJ^ ^njGjs 


(D12) 


For the kurtosis we need the fourth //-derivative. It can 
be straightforwardly derived from (D12) with (D9) and 
(iDlip, and leads us to 




dji"^ 

-f- \2,^^2iGif- ^^jGji T ^^jGji D^2^j 

^fijGjj- D^2^Gjs ^krs 
'^^fiiGik ^filGij-^fj^jGjs ^j_ikrs 

—190 o o o 

^ ^fll2 '^1232 ^ ^ ^k-3l32 

\ 2 ,^^i^Gi^i^^fj,i^iGik ^fijGji^fij^Gj^j^ ^kij 2 
311 ^ 21^2122 ^/i 23 1^23/4 11/224/G/fc ^^jGji ll/i_ 7 i^jij 2 ^kij 2 
T 11/222 Guk n^iG ir ^filGis 11/2j Gpn 

krst ■ (D13) 


We have also used (D8) and (D13) to compute higher mo¬ 


ments. A poten tial problem arises when the determinant 
of D(2) in ( |D7[ ) is vanishing, which hampers the inverse 
of and results in singularity for the propagator Gij. 
Therefore, computations based on this analytical method 
need further investigations. 
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FIG. 18. Quadratic (left panel) and quartic (right panel) baryon number fluctuations as functions of T calculated in the 
simplified PQM model, with the dependence of Polyakov loop on T as input. Here we compare two different truncations: LPA 
and full. 




T [MeV] T [MeV] 

FIG. 19. Left panel: kurtosis of the baryon number distribution in the simplified PQM model. Right panel: the red-dashed 
curve is just the same one in the left panel, in comparison to the blue-solid line, which is obtained with the Polyakov loop set 
to be zero artificially below T = 200 MeV, while unchanged above this temperature. See text for more details. 


Appendix E: Kurtosis in a simplified PQM model 

The representation of the higher moments in Ap¬ 
pendix suggests the potential importance of a self- 
consistent computation of all quantities, in particular 
that of the higher L, L-moments of the glue potential. 
Here we check this hypothesis in the simplified PQM 
model: the T-dependent Polyakov loop expectation value 
is used as an input, and its dependence on the chem¬ 
ical potential is neglected. As external input we use 
the Polyakov loop expectation value as computed in [52] . 
Note that the LPA computation with a full effective me¬ 
son potential there is very close to the present LPA com¬ 
putation within a Taylor expansion. 

Fig. |18| shows the quadratic and quartic baryon num¬ 
ber fluctuations in the simplihed PQM model within LPA 
and full approximation. One observes that the fluctu¬ 


ations in the improved approximation are larger, and 
vary more rapidly with the change of temperature dur¬ 
ing the chiral crossover, similarily to the results in the 
self-consistent computation, see Fig. [T3| 


The kurtosis in the simplified PQM model, however, 
looks similar to that in the QM model: the missing self- 
consistency of the glue background leads to a qualitative 
failure in the hadronic phase, see Fig.[^in comparison to 
the self-consistent results in Fig. M Instead, the kurtosis 
in Fig. is similar to that in the QM model, in agree¬ 
ment to an earlier mean field analysis in m, see Fig. 5 
and Fig. 7 there. If we enforce the hadronic nature with 
LL = 0 for T < Tc and LL = 1 for T > Tc the asymp¬ 
totic low and high temperature behaviour is restored, see 
the right panel of Fig. Note however, that this limit 
does not hold in the hadronic phase, instead we have (46) 
and (47). 
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